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We have performed inelastic neutron scattering on the near ideal spin-ladder compound 
La4SrioCu24C>4i , which will serve as a model system for the investigation of doped ladders. A 
key feature in the analysis was the separation of one-triplon and two-triplon scattering contribu- 
tions due to an antiphase rung modulation. This enabled extraction of the one-triplon dispersion 
which has a spin gap of 26.4HL3 meV, and for the first time two-triplon scattering is clearly observed. 



The exchange constants are determined by fitting the data to be J; e3 = 186meV and J r 



= 124 meV 



along the leg and the rung respectively. In addition to the surprisingly large difference between Ji eg 
and Jrung a substantial ring exchange of J cyc —31 meV was confirmed. 

PACS numbers: 78.70.Nx 75.30.Ds 75.10.Jm 75.30.Et 
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The two-leg ladder, consisting of two spin- 1/2 chains 
coupled by a Heisenberg interaction, can be thought of as 
a strip, two sites wide, cut out of a 2D square lattice anti- 
ferromagnet. The immediate effect of coupling the chains 
is to confine pairs of spin-1/2 spinons, and to generate 
a gapped spin singlet ground state. Short ranged singlet 
spin pairings on the rungs l/y/2{\ U) - I IT)} dominate 
the ground state, and in an inelastic neutron scattering 
measurement a triplet of well-defined spin-1 magnons, 
known as triplon [lj excitations, rather than a spinon 
pair continuum will appear. Further, strong higher-order 
ground state fluctuations should be manifest as multi- 
triplon continua. The spin pairings are particularly im- 
portant as they provide a mechanism for the formation 
of charged hole pairs when doped, giving rise to con- 
densed charge density wave or superconducting states 
Hi, SlU. Cuprate spin ladders share many features of 
the 2D high-2}, materials (pseudo-gap, non- Fermi liquid 
behavior) |7j, |8|, |9j , and have similar coupling parameters. 
Their excitation spectra extend to several hundred meV 
and the development of time-of-flight spectroscopy at 
spallation neutron sources allows access to the full spec- 
tra resolved in wave- vector and energy for the first time. 
Motivated by this, we have measured the one- and two- 
triplon spectra of a cuprate ladder in a Mott-Hubbard 
insulating (MHI) phase - La4SrioCu2404i- using inelas- 
tic neutron scattering (INS). A full description of the 
dynamics is achieved using the continuous unitary trans- 
formation (CUT) method. Electrostatic environments of 
the oxygen tunneling result in strongly dissimilar rung 
and leg exchanges, and this combined with a substantial 
cyclic charge fluctuations render the magnetic state of 
the system near quantum critical (gapless). 

The geometry of cuprate ladders in La4SrioCu2404i 
is depicted in the upper part of Fig. 1. Strong anti- 



ferromagnetic super-exchange occurs through the 180° 
Cu-O-Cu bonds with very weak exchange between lad- 
ders via the 90° bonds. An approximate electronic 
Hamiltonian for the cuprates is the one-band Hubbard 
model 

H = - ko ( c i^ c J> + H - c -) + U Y1 n '^ ni >l w 

with on-site Coulomb energy £/~3.5 eV, inter-site hop- 
ping between nearest neighbor pairs of t ~ 0.3 eV, and 
k = t/U ~ j2 10]. Charge fluctuations such as cyclic 
hopping processes, (suppressed in the k — > limit), are 
reflected in modifications to the Heisenberg spin Hamil- 
tonian with the cyclic hopping contributing a four spin 
interaction term. In fact computations in the more real- 
istic 3-band Hubbard model 11] (which include 02p and 
Cu3d orbitals) indicate only short range Heisenberg and 
four spin interactions are important with diagonal cou- 
pling across the CuO square plaquette and further neigh- 
bor interactions negligible. The resulting spin Hamilto- 
nian is 



H — Jrung ^ ' &i,2 + Jleg ^ &i,T ' Si 
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Hcyc — Jcyc J]] [ (Sj, 1 • Sj+i t i ) (Sj.2 1 Sj+1,2 ) + 
plaquettcs 

(Sj,i • S^2)(Si+i,i • Si+i,2j — (Sj,i • Si+i,2)(Si+i,i • Sj^)], 

The anticipated dispersion relation for the spin-1 triplon 
states is plotted in the lower part of Fig. 1. For the lad- 
der the relative ratios of the gaps at wave-vectors along 
the leg direction Q c = 0, 0.25, 0.5 (in units of ir/2c) are 
modified substantially by the strength of J cyc - 
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wavevector along the leg direction Q (2n/c) 



FIG. 1: (Color online) Upper part: Arrangement of the lad- 
der planes in La4SrioCu2404i. Big spheres represent the cop- 
per ions, small spheres oxygen atoms. The dashed rectangle 
displays the unit cell in the (010) crystallographic plane for 
the a axis and half of the unit cell in the c axis. Lower 
part: S=l excitation spectrum. There is a gap from the 
S=0 ground state to the S—l triplon band (red solid curve). 
The two-triplon continuum states bounded by the two blue 
dashed curves can be constructed via the scattering of the two 
trip Ions illustrated by the arrows. A typical wave- vector cut 
(horizontal) and energy cut (vertical) are indicated. 

Also shown are the boundaries of the two-triplon con- 
tinuum. The intensity of the two-triplon spectrum, 
wave-vector and energy distribution is strongly depen- 
dent on the triplon interactions and exact composi- 
tion of the ground state. Importantly, a bound spin-1 
state is predicted below the two-triplon continuum. The 
four spin exchange term frustrates the formation of this 
bound mode lowering the binding energy and for coupling 
strengths of order x cyc — J cyc / ' J rung ~ 0.25 will disrupt 
binding completely. 

A remarkable property of the ladder geometry is that 
the sectors with even and odd triplon number do not mix 
because of different parity w.r.t. reflection about the cen- 
tcrlinc of the ladder. Hence even and odd contributions 
can be Fourier decomposed by the rung coupling, and 
the neutron scattering matrix elements can be expressed 
as a product of a neutron structure factor (NSF) depen- 
dent only on the wave- vector along the leg direction q c , 
multiplied by a rung wave-vector q a modulation. The 
one- and two-triplon spectra are in anti-phase with rung 
modulations (1 — cos(g.a)) and (l + cos^.a)) respectively 
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FIG. 2: (Color) (a) Cut along the rung wave-vector shows 
the intensity modulation Upper part: for energy of 30±5 meV 
and leg wave- vector of 0.45±0.05 in the rung direction for the 
one-triplon scattering. Lower part: for energy of 300±10 meV 
and leg wave- vector of 0.25 ±0.05 for two-triplon. Solid line 
described in text, (b) Typical wave- vector cut through one- 
triplon along Q c for energy of 130 ± 5 meV. The solid line 
is a Gaussian fit to extract peak position, (c) Energy cut 
for 0.48 < Qc < 0.52 to determine the energy gap. Light 
grey points are the corresponding cut for a hole-doped ladder, 
(d) Experimental points from neutron data. Theoretical one- 
triplon dispersion curves shown in solid red and black. Black 
curve includes interladder coupling. Theoretical lower bound- 
ary of one-triplon scattering shown in solid blue. Dashed lines 
show respective scattering without cyclic exchange. 



[12J , and this can be used to separate them in a neutron 
scattering experiment. 

Here we investigate La4SrioCu2404i: Although intrin- 
sically hole-doped, X-ray absorption spectroscopy reveals 
that in La x (Sr, Ca)i4_ x Cu2404i with , residual holes are 
located in the more electronegative chains leaving 
ladders undoped [fH Il5l | . No triplon-hole scattering is 
observed in the heat conductivity [Tisl . It7| . 

Large La4SrioCu2404i single crystals were prepared 
using the ' Travelling Solvent Floating Zone '-method at 

9 bar oxygen pressure. Inelastic neutron scattering mea- 
surements were performed using the MAPS spectrome- 
ter at ISIS, Rutherford Appleton Laboratory, U.K. Three 
crystals of total mass 23 g were co-aligned with resulting 
mosaic spread of 0.4° in the b-c plane and 5.5° on the a-c 
plane. The samples were mounted with the (Okl) recip- 
rocal lattice plane horizontal, and the c axis (leg direc- 
tion) perpendicular to the incident neutron wave-vector 
hi in a cryostat giving a measurement temperature of 

10 K. The wave- vectors are labeled Q c along the leg and 
Q a along the rung. The magnetic cross-section was nor- 
malized using incoherent nuclear scattering from a vana- 
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FIG. 3: (Color online) (a) Energy cut for two-triplon at 
Q c = 0.5±0.1. Solid line is two-triplon model. Intensities had 
to be increased by 15 % to match data. Light grey points are 
the corresponding cut for a hole-doped ladder. Thin dashed 
line is the effect of absent cyclic exchange for fixed J Ttm g- (b) 
Relationship of uionset and Jrung w.r.t. x cyc . For an onset en- 
ergy of 225 meV Jrung = 124 meV and x cyc = 0.25 are needed 
(see text). J cyc — x cyc Jrung- (c) Frustration of bound mode 
with increasing cyclic exchange, unbinding difference between 
Uonset and energy of bound mode. 



dium standard. Data were collected for incident energies 
Ei— 75.55, 363.7 and 606.2 meV. The energy runs were 
monochromated using a Fermi chopper at 300 Hz for the 
low energy run and at 600 Hz for the other runs. 

Contributions mainly from one-triplon or two-triplon 
scattering were separated using the antiphase rung mod- 
ulation. The measured modulation is shown in Fig. 2 (a) 
with the one-triplon modulation fitted to (1 — cos{q ■ a)) 
and the two-triplon modulation to (1 + cos(g -a)), with 
corrections for the anisotropic magnetic form factor of 
copper 18| applied to the fitted functions. One-triplon 
data used in subsequent analysis was taken at the max- 
ima between 1 < | Q a \ < 2, see in Fig. 2 (a). For the two- 
triplon data, two regions of maximum two-triplon inten- 
sity < | Q a \ < 0.5 and 2.65 < | Q a \ < 3.35 were added tak- 
ing magnetic form factor corrections fully into account. 

The non-magnetic background correction for the one- 
triplon signal was deduced via interpolation of the inten- 
sity found at energies above and below the one-triplon 
dispersion. The background subtracted data is displayed 
in Fig. 4 (a) . For the two-triplon scattering, the back- 
ground was deduced from the scattering observed at 
Q c = 0, 1, 2 where the magnetic signal is expected to 
be weak. The result is shown in Fig. 4 (c). 

The one-triplon dispersion was determined by extract- 
ing positions and intensities of measured peaks by least- 
squares fitting of Gaussians. A typical cut is shown in 
Fig. 2(b), this gives Q c at energy of 130±5meV indi- 
cated in Fig. 1 by the horizontal slice. These along with 



data points from different cuts are displayed in Fig. 2 (d). 
Heisenberg couplings alone are unable to account for 
the energy gaps of 197, 286, 26.4 meV, at wavevectors 
0, 0.25, 0.5 respectively. Our computations using the 
CUT method [12J reveal that a substantial four-spin com- 
ponent J cy c= 31meV and Heisenberg exchange strengths 
along the legs Ji eg = 186 meV, and rungs J rung = 124 meV 
that are very different 19| are required to explain this. 
Computing the entire dispersion using these values (solid 
red line) results in excellent agreement with the data [2(J . 
The lower boundary of the two-triplon scattering was ex- 
tracted in the same manner and the data points along 
with the computed lower boundary of the two-triplon 
continuum (blue curve) are plotted in Fig. 2 (d) . Also 
shown is the computed one-triplon (and two-triplon) dis- 
persion in the absence of cyclic exchange (dashed curves) . 
One can see the necessity of introduction of the cyclic ex- 
change to describe the data within its error bars. 

The CUT method is also used to compute the NSF 
which agrees very well with the measured scattering, see 
Fig. 4(b) and (d). A complete simulation of detector 
coverage and multicrystal mount is included with the 
detailed angular mosaics of each crystal, measured via 
neutron Laue, accounted for using a Monte Carlo inte- 
gration. The uneven intensity profile of the data with 
wavevector is due to a combination of detector coverage 
and crystalline mosaic. Fig. 2(c) shows an energy cut 
at Q c = 0.5 (vertical slice in Fig. 1, wavevector range 
0.48 < | Q c \ < 0.52). The computed NSF including instru- 
mental corrections is shown as the red line and matches 
the data. Note the fitted gap is 26.4 ± 0.3 meV which is 
lower than found in previous measurements pH |22| . 

The effect of cyclic exchange on the two triplon binding 
is illustrated in Fig. 3. The data points in Fig. 3 (a) show 
a two-triplon wave-vector cut of Q c — 0.5±0.1; the red 
solid line is the same cut performed through the model 
with the determined exchange parameters. Fig. 3 (b) il- 
lustrates the dependence of the continuum onset energy 
uionset (measured as 224 meV), with J rung determined 
to match the one triplon gaps at Q c = 0.25, 0.5, on 
cyclic exchange ratio x cyc = J cyc / Jrung- Keeping J rung 
fixed and assuming no cyclic exchange a bound mode 
12meV below ui onset would appear, shown in Fig. 3(c) 
u>bindin g = 12meV. Such a shift is indicated by the black 
dashed line in Fig. 3 (a) and would be detectable within 
our resolution limits confirming we are on the binding 
threshold. 

Figs. 4 (c) and (d) show that nearly all the S = 1 two- 
triplon scattering weight is concentrated near the lower 
continuum boundary, and this is reproduced very well 
by calculation. Van Hove singularities in the density of 
states of phonon assisted IR absorption measurements 
indicate that 5 = bound states survive as they are 
more stable in quantum antiferromagnets |l2t l23t |24| . 
The CUT method also predicts the relative one and two 
triplon neutron scattering intensity scales. The scale fac- 
tors applied to the model above are within 15% of the 
computations, i.e. within experimental accuracy. 
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FIG. 4: (Color) (a) Measured one triplon data with non- 
magnetic background subtracted. The white curve gives the 
theoretically calculated one triplon dispersion curve, (b) Cal- 
culated one-triplon scattering with corrections for instrumen- 
tal resolution and magnetic form factor, (c) Measured two 
triplon data with nonmagnetic background subtracted. The 
white curve gives the lower boundary of the two-triplon dis- 
persion curve, (d) Calculated two-triplon scattering. 

The large difference between rung and leg exchange im- 
plies stronger hopping along the legs ti eg > t rung despite 
a longer pathway {di eg = 3.96A, d runn = 3.84A). Contri- 
butions from hopping between 02p orbitals |25j depend 
acutely on electrostatic environment 26]. For ladders, 
Fig. 1, 02p hopping along the leg, if' is enhanced by two 



neighboring Cu 2+ ions compared to the rung, tf pp which 
has one Cu 2+ neighbor [2(j accounting for the stronger 
tieg. Further, admixture of Cu4s states has been pro- 
posed to enhance i; eg hopping (by 35 %) [53, H3| com- 
pared to t rung . Comparison of our values of Ji eg , J run g, 
and J cyc with theoretical calculation should lead to a 
more refined model for cuprate electronic structure. 

Approximate parameters for (1), i; eg = 0.42, t rung — 
0.34, U = 3.72 eV, follow from low order perturba- 
tion theory [2^, but recourse to 3-band computations 
is necessary for high accuracy. To account for the 
CDW states shown at low temperatures in hole-doped 
ladders Sri4_ 2; Ca x Cu2404i |3G| electrostatic interactions 
Y,i, a =L,R Vleg n i,a.n>i+i,a + V rung rii, L n hR have been intro- 
duced in an extended Hubbard model [3l|]. Our MAPS 
measurements on the x = 11.5 doped ladder (p cx T above 
Tcdw ~ 80K) which has a modest doping of 5% holes on 
the ladder, show the one-triplon spectra to split into mul- 
tiple branches and for the gap to develop into a pseudo- 
gap (see 2(c)). Further the boundaries of the two-triplon 
spectra undergo dramatic changes and the lower bound- 
ary at Q c = 0.5 collapses (see Fig. 3(a)). Understanding 
the dramatic evolution of the spin correlations from the 
MHI to doped metallic ladder is a challenging problem, 
but is necessary to bridge the gap between understand- 
ing INS from insulating quantum magnets versus strongly 
correlated metals such as high-T c . 

In conclusion, our high resolution INS measurements 
establish the full one- and two-triplon spectra in the MHI 
cuprate ladder, and also the accuracy of the continu- 
ous unitary transformation method. Computation of the 
magnetic spectra in doped ladders within the framework 
of an extended Hubbard model is urgently needed. 
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